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Abstract 



Alternative definitions of the Born approximation and the distorted-wave 
Born approximation within the framework of the configuration-space Faddeev 
equations are explored. The most natural definition does not correspond to the 
Born approximation derived from the Schrodinger equation, even though the 
exact T-matrices for both formalisms are equivalent. The Schrodinger form is 
optimal, although it is shown that the differences are numerically unimportant. 
The DWBA corresponding to the Faddeev equations is not channel symmetric, 
although numerically this is unimportant for the p-d (Coulomb) case. Conver- 
gence of the Born approximation partial-wave series is briefly investigated for 
p-d and n-d scattering below breakup threshold. 
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1 Introduction 



Modern few-nucleon calculations are in principle a rich source of information about 
the nuclear force, provided that numerical results can be obtained with sufficient 
precision. The required level of precision depends sensitively on the observable being 
calculated. That level ranges from modest in the case of unpolarized observables such 
as the differential cross section to high in the case of some polarization observables, 
since the latter often depend on delicate cancellations between partial waves. One 
of the strengths of the few-nucleon field is that accurate calculations for realistic 
potentials are now routinefl, 2, 3]. 

One of the consequences of this numerical precision requirement is that a substan- 
tial number of partial waves may be required. This is not such a dire problem for 
momentum-space techniques, many of which generate a Born series for the T-matrix 
as the preferred method of solution. Such a series should converge faster as the angu- 
lar momentum increases. That is not the case, however, in configuration space, where 
each partial wave is solved in its entirety. Because more and more partial waves of the 
potential contribute as the angular momentum increases, this significantly increasing 
effort is rewarded by decreasing return as the size of higher-partial-wave T-matrix 
elements becomes smaller and smaller. One way to increase efficiency in the process 
is to use the Born approximation for higher partial waves, which significantly reduces 
the effort [4]. 

We detail below several pitfalls that arise in defining the Born approximation 
using configuration-space Faddeev techniques [5]. While none of these difficulties are 
very important in practical calculations (although this is not a priori obvious), it 
is nevertheless worthwhile to avoid them altogether. These problems are derived in 
Section 2 for the Born approximation (BA) and in Section 3 for the distorted-wave 
Born approximation (DWBA), and a few numerical results are given and discussed in 
Section 4, which is followed by our summary and conclusions. We restrict ourselves to 
energies below the deuteron-breakup threshold, and we ignore three-nucleon forces, 
which are an inessential complication. 



2 Faddeev Equations and the Born Approximation 



The Faddeev decomposition in configuration space [6] can be pictured most easily as 
dividing the wave function into parts corresponding to different asymptotic processes. 
For three (identical) particles this leads to three equations 
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where i,j,k take on the values of the cyclic permutations of 1, 2, and 3, and V* = 
V(xj). The Schrodinger wave function is given by 

* = ^i(xi,yi) +^2(x 2 ,y 2 ) +^ 3 (x 3 ,y 3 ) , (2) 

with arguments fashioned from the coordinates, r,, of all % particles: 

Xj = Tj - r k , (3a) 

and ^ 

y t = -(Tj + T k ) -Ti. (36) 

The coordinate Xj is the distance between particles j and k, while y» is the distance 
between particle % and the center of cluster j + k. Depending on the desired application 
a normalization factor of l/y/3 may be used in the definition of ^[7]. With these 
definitions the sum of the three independent Eqns. (1) yields the original Schrodinger 
equation 

(E -T -V 1 -V 2 -V 3 )^ = 0. (4) 

The coordinates (3) used in the Faddeev wave functions in Eqn. (4) generate different 
asymptotic configurations (or channels) that facilitate the appropriate boundary con- 
ditions for this partial-differential equation. Used in this way the Faddeev equations 
are merely another (equivalent) form of the original Schrodinger equation. 

The appropriate boundary conditions for elastic scattering (below breakup thresh- 
old) are easily seen from Eqn. (1). For large |yj| the scattered wave function factorizes 
into the product of a deuteron wave function 0d( x i) and an outgoing wave, x{Yi)i 
which guarantees that the right-hand side of Eqn. (1) vanishes in this limit. This 
product of 0d and a plane wave (or a Coulomb wave in case of p-d scattering), x(y»)> 
defines the "free" channel wave function, <^[7]: 

(f>i(xi, yi) = M*i)x(yi) , (5a) 

satisfying 

(E — T — Vi)(j>i = , (56) 
which we can use to construct the exact T-matrix and its Born approximation. 

The Faddeev T-matrix is most easily developed by projecting Eqn. (1) with 
The resulting left-hand-side integral vanishes because of Eqn. (5b), except for surface 
terms that constitute the T-matrix (just like the two-body case [8]) 

t F = (0ilW 2 + ^3>, (6) 

and any other permutation of Eqn. (1) yields the same result. This can be converted [9] 
to a more recognizable form by rearranging Eqn. (5b) to (0i|Vi = (4>i\(E — T), 
making use of the Hermiticity of the Hamiltonian for these wave functions below 
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breakup threshold, and employing Eqn. (1) to produce the Schrodinger form of the 
T-matrix[10] 

t s = (fa\V 2 + V 3 \^). (7) 

The potential V 2 + V 3 is the usual "inter-cluster" potential[10] between the nucleon 
projectile and the deuteron target that generates the scattering T-matrix. The two 
forms of the T-matrix given in Eqns. (6) and (7) are exactly equivalent if the wave 
functions are exact. 

Born approximations can be obtained from Eqns. (6) and (7) by replacing the ipi 

by (f)i and * by $ = 0i + 2 + 3 : 

t* A = (0iN0 2 + 3 >, (8a) 

and 

t| A = (<f>i\V 2 + V 3 \$). (86) 

In order to facilitate the comparison of the two Born approximations in Eqns. (8) it 
is convenient to introduce the compact notation 

(ijk) = (HVjlfa) . (9) 

Interchanging the labels of pairs of variables leads to (132) = (123), etc., while the 

trick used to convert Eqn. (6) to Eqn. (7) leads to (122)=(112), etc. We thus arrive 
at 

4 A = 2(112) , (10a) 

and 

tf A = 2(112) + 2(121) + 2(123), (106) 

which are distinctly different. One can also use these tricks to demonstrate that the 
Faddeev form of the Born-approximation T-matrix is channel symmetric: 

t f = + 03 1 ) = (02 + 03l Wl) = ^ , (IOC) 

where we have attached channel labels that were previously suppressed. In other 
words the "post" and "prior" forms of the Born approximation are identical, a result 
long known for £g[10]. The numerical significance of the difference of Eqns. (10a) and 
(10b) will be discussed in Section (4). 

An important point concerning the Born approximation has been made by Kievsky 
et al[4\. The Kohn variational functional can be written as [7] 

I(V) = (V\H-E\V) = (ij; 1 \H-E\ifr), (11) 

where the last step follows by interchange of variables and by removing the normaliza- 
tion factor l/\/3 from each wave function. Varying \I> (and -0i) about the plane- wave 
limit leads to 

tKohn = (0l|V2 + V 3 \$ ) + 0{V 2 ) , (12) 
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since any change in \l> beyond the plane-wave part must be proportional to the po- 
tential, and the Kohn error is proportional to the square of that change in \I>. For 
this reason tg A in Eqn. (8b) gives the optimal Born approximation. The two terms 
missing in tp A are "hidden" in ip2 +^3. Even though t^ A is linear in the potential V\, 
that potential serves partly to drive the dynamics (since Eqn. (6) is exact) and partly 
to bind the deuteron; the latter terms do not directly contribute to the T-matrix, 
and serve to "hide" V 2 and V 3 terms residing in fa and Thus the Faddeev Born 
approximation has first-order (in V) errors, unlike the Schrodinger result. 



3 Faddeev Equations and the DWBA 

The remaining task is to compare distorted-wave Born approximations from the Fad- 
deev and the Schrodinger equations. This is complicated by the fact that long-range 
forces (such as Coulomb) have been treated in the Faddeev formalism in a variety 
of ways, reflecting the non-unique character of the Faddeev equations themselves. A 
powerful formal device is the introduction of Faddeev distortion potentials, JQ. We 
can replace Eqn. (1) by [7] 

[E — T — Vi — Xj — X k }^ = (V, - + i> k ) . (13) 

If one sums the three independent (cyclic-permutation) members of this set, one 
obtains the original Schrodinger equation in Eqn. (4) (i.e., the Xi cancel), and \I/ is 
free of the arbitrary X it as well, although individual tpi are not free. It is effective 
to choose the X± so that long-range distortions are built into the wave functions ipi. 
As a specific example we treat the important case of long-range Coulomb distortions, 
and write Vi = V? 1 + Vp, where and Vf are the strong-interaction and Coulomb 
parts of the nuclear force. In this case we have found it most efficient numerically to 
choose Xi = V i , which places all of the Coulomb interaction on the left-hand side of 
Eqn. (13). Other choices are possible and have been used by others[ll, 12]. 

We begin our treatment of DWBA by noting that if a long-range force modifies the 
asymptotic form of the inter-cluster scattering wave function, x(y»)> our basis wave 
functions (i.e., the "free" scattering wave functions) are modified and correspond to 

(E-T-V 1 -U 1 )<f )l = 0, (14) 

where U\ is an appropriate long-range interaction potential (e.g., Za/y for p-d Coulomb 
scattering). We next rearrange Eqn. (13) (choosing % — 1) by moving all X-terms to 
the right-hand side, and then adding and subtracting a potential Lq: 

[E-T-V x - = (X 2 + X 3 - + (Vi - Xjfa + • (15) 

Projecting this equation on the left with gives a left-hand-side that vanishes 
because of Eqn. (14) except for surface terms (that comprise the T-matrix, as before) 

h = (0i |*2 + X 3 - Wi) + (0i I Vi - X x \i> 2 + ^3) , (16) 
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where we ignore here (and hereafter) the T-matrix coming solely from the potential 
Ui (e.g., the point-Coulomb T-matrix), which is important but peripheral to our 
discussion. The barred-t indicates asymptotic states modified by U±. Choosing Xi = 
Vp in the Coulomb case leads to our preferred form of the Coulomb-distorted-wave 
Faddeev T-matrix, tp : 

*F = (011^1^2 +^3> + (0l|V C - Wi) , (17) 

where V° is V 2 C + V 3 C , and we have used the fact that Vp vanishes. If we again use 
{4>i\Vi = (<f>i\(E — T) in Eqn. (16), we obtain after some algebra 

*s = i^ + Vs-U^), (18) 

a remarkably different-looking, but equivalent, result. Equation (18) is the usual 
Schrodinger form of the distorted-wave T-matrix, and differs from Eqn. (7) only by 
the Ui term [10]. 

Although the two forms of the T-matrix in Eqns. (16) and (18) are formally 
identical, they lead to very different DWBAs. Replacing ipi by 0, and ^ by $ in 
Eqn. (16) leads to a formally unacceptable DWBA, because it would depend on the 
unphysical (i.e., artificial) Xi 

f DWBA = (0i | X2 + X;j _ + _ Xi | 02 + ? (1Q) 

although the choice of Xi leading to Eqn. (17) obscures this fact. The DWBA 
obtained from Eqn. (18) 

i£ WBA = (<Pi\V 2 + V 3 - U^) , (20) 

is free of that disease because Eqn. (18) depends only on \1>, and not the individual 
ipi. The Schrodinger DWBA is the optimal one in the sense that corrections to it are 
second order in the potential. The techniques that led to Eqn. (10c) also demonstrate 
that the distorted-wave T-matrix in Eqn. (19) is not symmetric under interchange of 
incoming and outgoing channels (post and prior forms), because of terms involving 
X\ and U\. This is a serious disease, since it violates a symmetry, although we shall 
see in the next section that the effect is sufficiently small that it is not a significant 
practical problem. Equation (20) can be shown by the same techniques to be channel 
symmetric. 

In summary the "Faddeev" form of the DWBA has three formally serious prob- 
lems: (1) it leads to a non-symmetric T-matrix; (2) it leads to manifestly non-unique 
results (viz., the X-terms); (3) it leads to errors linear in the potential. How impor- 
tant these flaws are in a numerical sense is discussed next. 
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4 Numerical Results 



We briefly report two numerical results on: (1) differences between the Faddeev and 
Schrodinger forms of the Born approximation, together with the magnitude of the 
lack of symmetry of the post and prior forms of the T-matrix; (2) the fractional 
errors in the difference of the exact Faddeev and Born approximation results for both 
p-d and n-d scattering below deuteron-breakup threshold. 

The form of the Born approximations given by Eqns. (10a) and (10b) suggests 
that numerical differences should be small. Below the threshold for deuteron breakup 
the region of configuration space where scattering can occur is very limited. The 
deuteron bound-state wave function has a finite extent, Rj, while the nuclear force has 
a significantly smaller range, Ry- For this reason the integrand of the configuration 
(123), for example, is important only in the very small region where the nuclear 
potential in variable X2 overlaps with a deuteron wave function in variable xi and 
another in variable x 3 . This occurs only for very small |yi|, where the wave function is 
suppressed for large values of the angular momentum. By far the largest configuration 
will be (112), since V\ overlaps easily with 0d( x i), and this result has a much larger 
overlap with the deuteron wave function in variable x 2 . Numerical results bear this 
out; the (112) matrix element completely dominates. This immediately suggests that 
the lack of symmetry between (off-diagonal) T-matrix elements is also numerically 
small. 

The importance of the Born amplitudes for higher partial waves (in J, the to- 
tal angular momentum of the channel states) is indicated in Table I, where both 
the differential cross section, da/dQ, and analyzing power, A y , have been computed 
using the AV18 potential [13]. These observables were computed in two ways: ex- 
act calculations through J = 11/2 (supplemented by Born approximation results for 
higher waves where needed) and various comparison calculations described below. 
The maximum value of the magnitude of the fractional difference between the exact 
and comparison calculations was computed for the angular range 6 = 30° — 150°. The 
restricted angular range was used because the analyzing power gets small outside that 
range, and relatively large fractional errors in that regime would be meaningless. In 
the small table labelled Truncation the comparison calculation had its partial-wave 
series truncated at the labelled J (i.e., that J- wave and all higher ones were omitted), 
leading to some very large errors. If one adopts the criterion that observables should 
be calculated to 1% accuracy, one needs to keep all matrix elements through J = 9/2 
for the differential cross section, and through J = 15/2 for the analyzing power for 3 
MeV p-d scattering. 

In the tables labelled Born Substitution the comparison case involved substi- 
tuting the (Schrodinger) Born amplitude for the labelled wave and all higher ones 
through the maximum value of J = 15/2. The differential cross sections are accu- 
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Table I. Maximum absolute values of the fractional differences between exact and 
comparison calculations for differential cross sections and analyzing powers versus 
partial wave. The comparison labelled Truncation deletes the labelled partial wave 

and all higher ones. The comparison labelled Born Substitution substitutes the 
Born approximation for the exact T-matrix elements for the labelled partial wave 
and all higher ones. Results are presented in the format x.x[n] = x.xlO n . 



Partial Wave 


5 7 
2 2 


9 
2 


n 

2 


13 
2 


15 
2 




Truncation — 


3 MeV p-d 








A{da/dfl)/(da/dQ) 

A(Ay)/Ay 


4.9[-l] 1.1[-1] 
3.1[ 1] 4.4[ 0] 


2.8[-2] 
1.5[ 0] 


7.8 [-3] 
3.9[-l] 


2.0[-3] 
8.6[-2] 


3.9 [-4] 
1.8 [-2] 



Partial Wave 


5 
2 


7 
2 


9 
2 


n 

2 


Born Substitution 


- 3 MeV 


n-d 




A(da/dQ)/(da/dn) 

A(Ay)/Ay 


1.4[-1] 
4.1[ 0] 


2.2[-3] 
2.2[-l] 


1.2[-4] 
1.9 [-2] 


8.9[-6] 
1.9[-3] 


Born Substitution 


- 1 MeV 


p-d 




A{da/dtt)/{da/dtt) 

A(Ay)/Ay 


7.0 [-2] 
8.6[ 0] 


7.1 [-4] 
1.3[-1] 


1.3 [-5] 
1.1 [-2] 


4. 3 [-6] 
2.7[-3] 


Born Substitution 


- 3 MeV 


p-d 




A(da/dn)/(da/dtt) 

A{Ay)/Ay 


8.5[-2] 
4.2[ 0] 


2.1 [-3] 
1.4[-1] 


9.7[-5] 
1.3[-2] 


5.8[-6] 
l.l[-3] 



rately computed (using the 1% criterion) by keeping exact matrix elements through 
J = 5/2, while the analyzing power requires exact J = 9/2 matrix elements. Exact 
configuration-space Faddeev calculations for large values of J (such as 11/2) require 
significant computational resources, and should be avoided if possible. 

The difference between the Faddeev Born approximation and the Schrodinger Born 
approximation depends on the particular matrix element, but typically is only a few 
per cent, and this is not likely to be very important. The Schrodinger form of the 
BA and DWBA is nevertheless demonstrably closer to the exact matrix elements for 
higher partial waves than is the Faddeev form of those matrix elements. The channel 
asymmetry in p-d scattering matrix elements at 3 MeV is typically a few times 10 -3 , 
and hence is numerically unimportant for the Coulomb DWBA problem. 

In summary, we have shown that the most natural Born approximation for the 
Faddeev T-matrix differs from that of the Schrodinger equation. We have shown 
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that the Faddeev DWBA has three formally serious properties: (1) it leads to a non- 
symmetric T-matrix; (2) it leads to manifestly non-unique results; (3) it leads to 
errors linear in the potential. In practical (i.e., numerical) terms these defects are 
not very important. Our analysis demonstrates that the dominant part of the Born 
approximation is contained in both approaches, and they differ only in much smaller 
contributions. 
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